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We provide a comprehensive picture of the jamming phase diagram by connecting the ather- 
mal, granular ensemble of jammed states and the equilibrium fluid through the inherent structure 
paradigm for a system hard discs confined to a narrow channel. The J-line is shown to be divided 
into packings that are thermodynamically accessible from the equilibrium fluid and inaccessible 
packings. The J-point is found to occur at the transition between these two sets of packings and is 
located at the maximum the inherent structure distribution. A general thermodynamic argument 
suggests that the density of the states at the configurational entropy maximum represents a lower 
bound on the J-point density in hard sphere systems. Finally, we find that the granular and fluid 
systems only occupy the same set of inherent structures, under the same thermodynamic conditions, 
at two points, corresponding to zero and infinite pressures, where they sample the J-point states 
and the most dense packing respectively. 



The way particles pack together influences the thermo- 
dynamic and mechanical properties of a wide variety of 
materials, including crystals, amorphous glasses and liq- 
uids, as well as many athermal materials such as granular 
sand piles and foams [l[ . Thermal and athermal materi- 
als have also been shown to have similar dynamic proper- 
ties as the point at which the particles jam together into 
a mechanically stable packing is approached. The jam- 
ming phase diagram 0, which locates the jamming 
transition as a function of various thermodynamic and 
dynamic variables, was introduced in an attempt to de- 
scribe the jamming phenomena of these systems within 
a unified conceptual framework. However, there is still 
considerable debate concerning the details of the jam- 
ming phase diagram 0, and in particular, how the 
J-point is connected to the thermodynamics and dynam- 
ics of these systems. For example, the J-point has been 
connected with the random close packing (| , while mean 
field studies , along with others [llj , suggest there 

a is a range of jamming densities, associated with a J-linc, 
starting from the lowest density jammed packing related 
to the J-point and increasing up to the most dense, glass 
close packing state. 

The goal of this letter is to explore the key elements 
of the jamming phase diagram through the analysis of a 
simple system for which both the statistical mechanics of 
the granular system, the equilibrium fluid and the con- 
nection between these two states of matter can be treated 
exactly within the inherent structure paradigm [I3 - ll6j |. 
This allows us to move beyond the mean field results 
which are restricted to the study of jammed packing pro- 
duced by quenching fluid states from above a density, <fid, 
where the fluid phases splits into an exponential number 
of states and particles are unable to diffuse out of their 
respective cages. Our results show that the J-line extends 
over a broader range of densities than expected. Further- 



more, the J-line is divided into packings that are ther- 
modynamically accessible or inaccessible from the equi- 
librium fluid. We also find that the J-point occurs at the 
maximum of the configurational entropy and marks the 
transition between the accessible and inaccessible pack- 
ings. 

We begin by generating the full ensemble of jammed 
states for a system of N, two dimensional (2D) hard discs 
of diameter a, confined to a narrow channel by two hard 
walls separated by a distance H/a < 1 + \/3/4. In D- 
dimensions, a spherical particle is locally jammed if it 
has at least D + 1 rigid contacts arranged such that they 
are not all within the same hemisphere. However, the 
local jamming of all the particles in a structure is a nec- 
essary, but not sufficient, condition to ensure collective 
jamming because the concerted motion of a group of par- 
ticles may allow the structure to collapse 17|, |l8| . The 
confinement of the present model prevents the discs from 
passing each other, which eliminates the possibility of 
collective motions of particles unjamming the packings 
and allows us to count all the collectively jammed struc- 
tures by simply considering local packing constraints. It 
also ensures that each disc can only interact with its 
nearest neighbors and the wall. As a result, there are 
only four particle configurations that satisfy the local 
jamming constraints, two dense configurations (denoted 
1 and 3) and two open, defect type configurations, de- 
noted 2 and 4 (see Fig. Q}. The volume associated with 
each configuration is given by Hl^j, where the is the 
longitudinal distance between neighboring discs' centres, 
with Z M = Z Ji3 = [H(2a - H)] 1 / 2 and l i<2 = l lA = o. 
Any jammed configuration can now be identified by a 
list of neighboring bonds. However, it is important to 
note that configurations containing neighboring defects 
are not allowed as the central disc in the local arrange- 
ment is unjammed. Furthermore, all the jammed states 
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FIG. 1. Analytical quench connecting equilibrium fluid con- 
figurations of discs i, m, n, j to the most dense (1,3) and least 
dense (2,4) local packing arrangements. Equilibrium configu- 
ration of 4 discs are initially mapped to a tangent discs config- 
uration by compression in the x-axis. The local arrangement 
of the discs contained in the product of the kernel K(yi,yj) 
and the dot product of vectors (Eq [5j maps the central discs 
to their jammed configuration. See text for more details. 



of the model are isostatic [19j . 

Following the statistical mechanics app roach to gran- 



ular materials suggested by Edwards [19|, [2Cj , the quasi- 
one-dimensional nature of the model allows us to use 
the transfer matrix approach to construct the ensemble 
of jammed states, taking bond incompatibilities into ac- 
count [2l|. For a fixed N, the volume of the system will 
fluctuate depending on the number of type 2 and 4 states 
in the configuration so we introduce a longitudinal pres- 
sure Pl-, as a conjugate to the volume, and fix the system 
at a constant temperature, T. The transfer matrix then 
takes the form: 



M 



/0 M 1>2 M 1>3 \ 

M 2>3 

M 3 ,i ' M 3A 

\Ma,i M 4 , 3 J 



(1) 



where My = Cy exp(— /SP^/io^y ), where ho = H — a. 
The exponential term is the Gibbs measure appropri- 



ate for the N,Pl,T ensemble and Cy is zero when the 
two bonds are incompatible and one otherwise. (3 = 
1/kT, where k is the equivalent to the Boltzmann con- 
stant in the jammed ensemble. In the thermodynamic 
limit, the partition function for the system is given by 
A(N,P L ,T) = NkTln(X), where A is the largest eigen- 
value of M. The configurational entropy, which is di- 
rectly related to the number of jammed packings, Nj, by 
S c = k In Nj((f>j), and 4>j are given 



S c /Nk = In A + T<9(ln A) /dT . 



AHL,, 4kTd(\n\)/dP L 



(2) 



(3) 



Fig. [2] shows S c /Nk as a function of <fij for the case 
when H/a = 1.866, obtained from a parametric plot 
of Eqs. [3] and [5] with respect to P^. The same dis- 
tribution of states for this system can be obtained us- 
ing a combinatorial approach [22[ and gives S c /Nk = 
(1 -0)ln(l -9) - 9 In 9- (1 - 26»)ln(l - 20), where is 
the mole fraction of defects (type 2 and 4 bonds). The 
advantage of the present method is that it allows us to 
follow how the granular system explores the packing land- 
scape as a function of the externally applied pressure and 
we have plotted the full equation of state (EOS) in Fig [3] 
While there is no internal pressure from the particles in 
a granular system, it is still necessary for it to do work 
against Pl if the system expands to samples less dense 
states, but there are fewer high density basins, so the bal- 
ance between these two competing elements results in the 
"equilibrium" condition for the granular system. In the 
limit that Pl^o oo, S c /Nk — s- as the system moves 
towards the most dense state with c/>j — 0.842, while as 
Pl^o — > 0, the system samples the jammed states asso- 
ciated with the maximum in S c where = 1/2- Vb/10 
and 4>j = 0.659. If we allow the pressure to go below 
zero, we can sample less dense packings and the system 
enters the least dense jammed state with </>,/ = 0.561, 
S c /Nk = 0, as Piho — > — oo. However, the hard particles 
have no attractive component to their interaction poten- 
tial that could sustain a negative pressure, suggesting the 
packings below the S c maximum are thermodynamically 
inaccessible. 

We now examine how the equilibrium fluid for this sys- 
tem samples the landscape. The exact partition function 
for this model was originally solved by Barker [23J, but we 
make use of transfer matrix solution by Kofke et al [24[ ■ If 
the positions of the discs are fixed in the y-direction, the 
configurational integral in the x-direction can be simply 
treated as a Tonks gas [2^]. Using a Laplace transform, 
the volume dependence is removed, leaving the transfer 
integral for the partition function Z in the NPT ensem- 
ble: 
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FIG. 2. S c /Affe versus <pj for ij/a = 1.866 showing the J- line. 
The cyan line highlights the set of <f>j belonging to the ac- 
cessible jammed states which is bracketed by the most dense 
(black point) and the J-point (green point), which also co- 
incides with the maximum in S c - The red line indicates the 
inaccessible packings leading to the least dense packing (red 
point). 



FIG. 3. Equations of state for H/a = 1.866. EOS for granular 
ensemble plotting pPhho as a function of 4>J (blue dot-dash 
line). EOS for equilibrium thermal fluid plotting /3Pl/io as 
a function of <j> (black solid line). <j>j of basins sampled by 
equilibrium fluid (red dashed line). EOS for fluid compressed 
using the LS method with different da/dt (dotted lines). (j>d ~ 
0.48 obtained from MD simulations. Insert: <j>j as a function 
of da/dt. 



Here A is the thermal wavelength, Pj, is the longitudinal 
pressure and K(y 1 ,y 2 ) = cxp[-P L h L x (y 1 ,y 2 )}, with y x 
and y 2 being the y-coordinates of two adjacent discs in 
contact. L x is the projection of the distance between the 
two contacting discs along the x-axis and is a function 
of 2/i 1 2/2 ■ Solving the eigenvalue problem associated with 
Eq. 2] [H[ yields the equilibrium equation of state for the 
fluid plotted in Fig. H 

To quench the equilibrium fluid to its inherent struc- 
ture, we take advantage of the information contained 
within the matrix K regarding the geometry of adjacent 
tangent discs. Starting from an equilibrium configura- 
tion, we translate the discs along the x-axis only so that 
discs are in contact with their nearest neighbors. Fig. Q] 
shows that the type of bond (1,3 or 2,4) between the two 
central discs (mn) that will result from further compres- 
sion can be determined from the sign of the product of 
areas made from the triangles created by particles i,m,n 
and m, n, j. The geometry of the four discs is contained in 
the chain product matrix K(y i ,y m )K(y m ,y n )K(y n ,y : j). 
The product area-vector-product rule, for triangles Aj m „ 
and A mn j, that determine the nature of the bond is, 

A 

imn • A nin j > bond mn irnri) — Ik.i 
Ai mn • A mn j < bond mn A?j? (mn) = a. (5) 

We can now define a new transfer matrix G, whose 
elements are weighted by the bonds they would jam to 



under the jamming criterion, 

i—ii ,z 4 

(6) 

For a system with periodic boundary conditions and A^— 2 
particles, the volume is given by, 

V™\= hm dLog[Tr(G)]/d 7 . (7) 

7->0 

This method represents an ideal, infinitely- fast, non- 
equilibrium quench of the liquid to its inherent structure 
and amounts to a Stillinger-map, where every configu- 
ration is mapped uniquely to the jammed state in the 
basin. Consequently, each jammed state is weighted by 
the volume of configuration space in the basin. Fig. [3] 
shows both the equilibrium EOS for the fluid as a func- 
tion of (f> and the <f>j of the inherent structure sampled by 
the fluid. We can now compare how the equilibrium fluid 
and jammed system sample the same inherent structure 
landscape as a function of the externally applied pres- 
sure. While they exhibit the same general trend by mov- 
ing to deeper inherent structures as Pl increases, the two 
systems visit different basins at the same pressure. The 
equilibrium properties of the fluid arise from competi- 
tion between the configurational entropy (the number of 
basins) and the free volume of a single basin, while the 
jammed ensemble has no free volume contribution. As 
a result, the fluid samples lower density basins than the 
jammed system at the same pressure. 
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However, the fluid and the jammed ensemble do sample 
the same set of basins at two particular values of Pl ■ In 
the limit that Pl — > oo, they both move towards the most 
dense state. There is only a single most dense structure, 
so the configurational entropy goes to zero, but there is 
no ideal glass transition in the fluid as this only occurs 
at Pl — > oo. The two systems also sample the basins 
associated with the maximum in S c as Pl — > 0, i.e. as 
the fluid approaches the ideal gas state. 

A simple thermodynamic argument can help us un- 
derstand which jammed states are thermodynamically 
accessible in hard sphere systems in general. The en- 
tropy of a hard particle system at </>, can be written as 
S f ((t>)/Nk = lnN g (ct>j) +]nQ g (<f>,cl>j), where Q g (<f>^j) 
is the configuration space of a single inherent structure 
basin that maps to <pj. At equilibrium, the system sam- 
ples the set of basins with <fij that satisfy, 

(dS f /d<f>j)t = . (8) 

The free volume equation of state 1261] a nd the equa- 
tion of state for constrained glasses [la, [iff suggest that 
(dQ g (</>,4>j)/d<f>j)4 > > 0. On the other hand, we gen- 
erally expect there to be a single maximum in N g (<f)j), 
with only a few very high or very low density states, 
because of the coupling between density and structural 
order [27| and estimates of hard sphere packing distribu- 
tions generally yield Gaussian type distributions jljj Gil ■ 
Combined, these conditions necessarily imply the equi- 
librium system is unable to sample basins with a <f>j 
lower than those at the S c maximum. If the system 
did sample states with <pj lower than the maximum, 
we would have (dN g (<j)j)/d4>j) > and the equilib- 
rium condition could never be satisfied. If the ideal gas 
does sample the associated with the S c maximum, then 
(dQ g (<f),<pj)/d(pj)^ = 0, which suggests that the basins, 
for all the glasses become the same in the low density 
limit. This is true for the quasi-one-dimensional system 
studied here, and for the one-dimensional system of non- 
additive hard rods, where at <$> = every particle in the 
every glass becomes a point cag ed by its two neighboring 
points on the line. Speedy |28| also found, for a binary 
mixture of hard discs in two dimensions, that the differ- 
ence in entropy between the free ideal gas and the ideal 
gas constrained to a single basin is independent of the 
basin's (j>j. However, a recent study [29| of disc pack- 
ings involving a small number of particles suggest the 
the basin volumes of individual packings have complex 
geometries at low densities suggesting it is not obvious 
that all packings will have the same basin volume in the 
ideal gas state. Nevertheless, Eq[8j in the limit 0^0, 
defines the lowest density jammed packing accessible to 
the equilibrium fluid, which also provides a thermody- 
namic definition of the J-point. 

We now compare our ideal inherent structure mapping 
with a more traditional compression scheme applied to 



hard sphere particles. We performed molecular dynam- 
ics simulations of our system using N = 10000 discs. 
Starting from a random, low density configuration with 
<fi = 0.05, the discs were compressed using a modified ver- 
sion of the Lubachevsky and Stillinger [30( scheme that 
ensures H/a remains constant as the discs are expanded 
at a rate of da/dt. The EOS for the system under differ- 
ent compression rates are plotted in Fig [3] The compres- 
sions follow the equilibrium EOS at low densities because 
the compression scheme allows the system to move be- 
tween basins as it evolves. Eventually the fluid falls out of 
equilibrium at higher densities when caging effects start 
to become important and the system becomes trapped in 
a glassy state consisting of a single basin on the inherent 
structure landscape. Continued compression leads to a 
jammed state where the pressure diverges. The jamming 
density of the glasses as a function of compression rate 
(see insert Fig[3]) was obtained by counting the number of 
defects in the glass and taking averages over 20 indepen- 
dent runs at each da/dt. Not surprisingly, the jamming 
density increases with decreasing da/dt as slower com- 
pression rates allow the system to remain in equilibrium 
longer. To reach the most dense state, the system would 
have to be compressed sufficiently slowly that it was in 
equilibrium at all points along the trajectory. Finally we 
note that the density at which the system starts to fall 
out of equilibrium as a result of compression in Fig.[3]rep- 
resents the fid used in the mean field approaches. The 
(f>js of packings generated from fluid configurations at 
this density are approximately 0.72, but our exact results 
show that there are many more jammed configurations 
with a lower <f>j. 

To conclude, our results show that the behavior of 
fluids, glasses and jamming can all be related through 
the underlying inherent structure landscape. Our model 
gives rise to a distribution of packings that is essentially 
the same as that speculated to exist by Ciamarra et al [9( . 
However, we also find that the fluid and jammed ensem- 
bles are connected thermodynamically on the jamming 
phase diagram at two points corresponding the least and 
most dense states accessible to the system from the equi- 
librium fluid. The least dense states accessible to the 
system occur at Pl = and are associated with the 
maximum in S c . These states are connected to ideal gas 
states by an ideal, infinitely fast compression and mark 
the transition between the thermodynamically accessible 
and inaccessible packings. The most dense state occurs 
at Pl = oo, with S c = 0, and is only accessible from 
an ideal gas starting configuration by an infinitely slow 
compression. 
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